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We present a comprehensive study of heat transport through small superconducting point con- 
tacts. The heat current for a temperature biased weak link is computed as a function of temperature 
and barrier transparency of the junction. The transport of thermal energy is controlled by the quasi- 
particle transmission probability for the point contact that couples the superconducting leads. We 
derive this transmission probability and results for the heat current starting from nonequilibrium 
transport equations and interface boundary conditions for the Keldysh propagators in quasiclassical 
approximation. We discuss the thermal conductance for both clean and dirty superconducting leads, 
as well as aspects of the nonlinear current response. We show that the transmission probability for 
continuum quasiparticle states is both energy- and phase-dependent, and controlled by an interface 
Andreev bound state below the continuum. For high transparency barriers the formation of a low- 
energy bound state, when the phase is tuned to (f> = n, leads to a reduction of the heat current 
relative to that for 4> = For low-transparency barriers, a shallow Andreev bound state just below 
the continuum edge is connected with resonant transmission of quasiparticles for energies just above 
the gap edge, and leads to enhanced heat conductance as the temperature is lowered below the 
superconducting transition. 

PACS numbers: 74.25.Fy,74.40.+k,74.45.+c,74.50.-|-r 



I. INTRODUCTION 

In recent years theoretical and experimental investiga- 
tion of quantum transport through small superconduct- 
ing hetero-str uctu res has increased significantly (see for 
example Refs. Ill2l and references therein). These struc- 
tures vary from length scales of the order of a microm- 
eter down to the atomic scale. Experimental progress 
is driven in part by the development of atomic-scale 
break junctions, sub-micron- to nano-scale lithography 
and the fabrication of novel multi-terminal mesoscopic 
structures, as well as experimental techniques and meth- 
ods of measuring the local temperatures in small meso- 
scopic structures 1^ 

Considerable effort has gone into the study of charge 
transport in mesoscopic devices, both the non-dissipative 
superconducting response as well as the dissipative 
current response in voltage-biased junctions^ Much 
less is known about the electronic heat transport in 
tempcraturc-biascd junctions. Here wc provide a detailed 
report of theoretical results on thermal transport in small 
temperature- and phase-biased Josephson junctions and 
weak- links as a function of barrier transparency, D, tem- 
perature, T, and phase bias, 0. We provide a deriva- 
tion of the transport properties of superconducting point 
contacts starting from the Keldysh formulation of the 
nonequilibrium transport equations for the quasiclassical 
Green functions, and extend our brief report in Ref. 
of the thermal conductance (linear response) to the non- 
linear current response. We use the Ricatti formulation 
of the tr ansp ort equations and follow closely the notation 
in Refs. |3i 

Small Josephson junctions are ideal systems for in- 
vestigating quantum effects on transport under well de- 
fined non-equilibrium situations. Here we focus on non- 



equilibrium transport induced by a temperature bias AT 
across the junction. The thermal current is shown to 
depend on the phase difference, 4>, across the junction. 
To quasiclassical accuracy, i.e. to leading order in the 
small parameters s = {Tc/ef,l/kf^o,...), the thermo- 
powcr is generally negligible, and in this limit there 
is no temperature-induced voltage across the junction. 
The phase remains constant, and thus the temperature- 
biased junction is in a stationary state. But, in gen- 
eral the thermo-power is non-zero, so to be more pre- 
cise we estimate the relative magnitude of the dissipa- 
tive charge current to thermal current induced by the 
temperature bias to be of order ~ s(e/fcs) where 
s ^ (Tc/ef) <C 1 is the relevant particle-hole asymme- 
try parameter. In a closed circuit particle-hole asym- 
metry leads to a temperature-induced voltage of order 
eV ^ kB^T{Tc/ef). Thus, for measurements of the 
heat current on a timescale t t ^ (tt / AT) {e f /Tc) , 
the phase is constant and we may consider the ther- 
mal current response in a stationary situation. This 
is in contrast to the dissipative charge current for a 
voltage-biased junction, where the voltage is generally 
so large that the a.c. Josephson effect develops. The 
charge current is then composed of a phase-independent 
d.c. current and the a.c. Josephson current which os- 
cillates with frequencies nwj, where ujj = 2eV/h is the 
Josephson frequency and n is an integer number liS. Thus, 
the temperature-biased junction opens up the possibil- 
ity of studying quasi-stationary, phase-dependent non- 
equilibrium transport processes. 

At a junction between two superconductors, there are 
in general bound states in the quantum well formed 
between the superconducting order parameters on the 
two sides^ In a clean system, without normal backscat- 
tering at the junction {D = 1), the formation of An- 
dreev bound states (ABS) can be understood in sim- 
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pie terms by considering a quasiclassical path: a closed 
path involving one electron-leg and one hole-leg, the two 
legs connected by Andreev reflections at the two sides. 
Since the superconducting phase is encoded in the co- 
herence factor during Andreev reflection, the resulting 
bound state energies depend on the phase. For con- 
tacts in which the distance between the electrodes is 
small compared to the coherence length, there is only 
one pair of bound states, with energies that disperse 



with the phase as ~ ±Ay 1 — D sin^(0/2), where 
< Z? < 1 is the transparency of the interface barrier. 
Since an electric charge 2e is transferred during Andreev 
reflection, the bound states participate in charge trans- 
port. In fact the phase dispersion determines the elec- 
tric supercurrent-phase dependence through the relation 
/ = (2e/h)deB/d4'- The continuum quasiparticle states 
do not contribute to the charge current; the bound states 
determine the current-phase characteristics. 

The situation is quite different when we consider heat 
transport. As is well known, Andreev reflcctioniii re- 
sults in strong suppression of the heat current through a 
normal-superconducting (NS) interface; an electron and 
a retro-reflected hole have equal, but oppositely directed, 
probability currents. Thus, only a fraction of the contin- 
uum quasiparticle states above the gap which are not 
retro-reflected contribute to the heat current across an 
NS interface. For two superconductors coupled by a 
point-contact wcaklink (ScS) the local excitation spec- 
trum near the contact plays a central role in regulat- 
ing the heat current through the interface. Although 
the continuum states above the gap carry the heat, the 
transmission probability of these excitations is strongly 
influenced by an interface Andreev bound state (ABS) at 
the point contact. Since the binding energy and spectral 
weight of the ABS are controlled by the relative phase 
of the two superconductors, the transmission probability 
and the resulting heat current are also phase dependent. 
For high transmission barriers the reduction of the con- 
tinuum density of states by the formation of the ABS for 
4> ^ leads to a suppression of the heat current compared 
to the case with = 0. For lower transmission barriers 
D < 0.4 and (/> 7^ 0, resonant transmission of quasiparti- 
cles with energies just above the gap leads to an increase 
in the heat transport relative to the normal state at Tc, 
over a wide temperature region below In this paper 
we investigate this resonance effect and the heat current 
in detail, including the dependence of the heat transport 
on the relative phase, the barrier transparency, the bar- 
rier model, impurity scattering, as well as the sensitivity 
of the heat transport to a finite temperature bias AT, i.e. 
beyond the linear response limit for the conductance. 

The paper is organized as follows. In Sec. ^ we de- 
scribe our model for small Joscphson junctions and point- 
contact weak links. In Sec. lIIII we derive the heat current 
for ballistic leads. We present the results for the linear re- 
sponse m Sec. lEI and the non-linear response in Sec.lvl 
We also discuss the tunnel limit D <C 1 in some detail, 
and examine the singularity that is encountered within 
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FIG. 1: Geometry of the temperature-biased Josephson con- 
tact of radius a. There is a potential barrier at 2 = with 
transparency D{pf), which connects trajectories with posi- 
tive and negative projection of the Fermi momentum on the 
z-axis. The labels indicate which Ricatti amplitudes are com- 
puted along each trajectory. The shaded boundaries define 
the region (~ ^ a) where superconductivity and the exci- 
tation spectrum are strongly modified for 4>i ^ 4>2- 



the tunnel Hamiltonian (tH) method that was previously 
used to calculate the heat conductance of an SIS tunnel 
junctioniiiiiSii^ We include a discussion of the heat cur- 
rent based on the tH method in the Appendix. Finally, 
the effect of disorder in the superconducting leads is ex- 
amined in Sec. IVII for superconductors in the diffusive 
limit. 



II. MODEL 

We consider two superconductors, denoted 1 and 2, or 
L (left) and R (right), connected through a small con- 
striction of diameter 2a at z = 0, see Fig. ^ The con- 
striction is assumed to be much smaller than the super- 
conductor coherence length. 



2a <^U = ^r 



(1) 



as well as the elastic mean free path, I, and inelastic 
mean free path, We will consider the ballistic limit, 
^ ^ ?A in Section IIIII and the diffusive limit ^ ^ 
in Section IVll In all cases, we assume that the inelastic 
mean free path is much larger than any other length scale 
in the problem. 

The potential barrier at z = (the vertical line in 
Fig. ^ is characterized by a transmission probability, 
D{\)f) for normal-state quasiparticles with Fermi mo- 
mentum Tpf incident on the interface. The reflection 
probabilty is given by R{pf) = 1 — -D(p/). The angle- 
dependence of D(p/) depends on the microscopic barrier 
and is not very important for our purposes; we use a 
(5-f unction barrier which gives. 
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Dp cos^ e 

1 - Do sin^ 



(2) 
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where 9 = arccos(v/-z) is the angle between quasiparticlc 
trajectory and interface normal, is the Fermi velocity 
and Dq is the transmission probability for quasiparticles 
incident normal to the interface. We develop the theory 
for heat conduction for arbitrary transparency of the con- 
striction: < Do < 1- This theory describes a variety 
of Joscphson weak links, including the tunnel junction 
{Dq ^ 1) and the pin-hole point contact {Dq — 1). 

For the order parameters of the two superconductors, 
we restrict ourselves to the singlet pairing case, but keep 
the orbital symmetry arbitrary. Thus, the order param- 
eter is in general assumed to be of the form Z(TyA(p^), 
where ay is the second Pauli matrix in spin space and 
A(p/) is the momentum dependent gap amplitude. The 
results reported below for the heat transport are for 
low- Tc superconductors with isotropic order parameters. 
However, the formalism is applicable to a broader range 
of Joscphson contacts, including anisotropic or unconven- 
tional spin-singlet pairing states. We note that current 
conservation is only guaranteed if the propagators and 
self-energies (order parameter, impurity self-energy, etc.) 
are calculated self-consistcntly. In the point-contact ge- 
ometry we can neglect the back-action of the current on 
the self-energies to lowest order in the small parameter 

The electrical and thermal resistances of the junction, 
both oc a~^, are much larger than the corresponding re- 
sistances of the leads. Thus, the phase change, as well 
as the temperature change, occur essentially at the junc- 
tion. We then take the local phases and temperatures 
near the contact to be equal to the reservoir phases and 
temperatures. Thus, we write the order parameters as 
Ai(Ti)e''^i and A2(T2)e^'^2 in superconductor 1 and 2, 
respectively. Here are the local temperatures near the 
contact given by Ti and T2, respectively. We define the 
phase difference over the junction </> = 02 ~ </>! and the 
temperature bias AT = T2 — Ti. In the linear response 
we will write Ti = T and T2 = T + AT. 



III. BALLISTIC CASE 

To study the transport of quasiparticles through su- 
perconducting point contacts as described above, we 
use the method of nonequilibrium quasiclassical Green 
functions li^iiSiiLiSiiS In this formalism the advanced 
and retarded Green functions, describe the lo- 

cal spectrum of excitations of the system, while the 
Keldysh Green functions, 25^, carry information about 
the nonequilibrium population of these states. The prop- 
agators ©^'^'^(p/, R; e, are 4x4 matrices in Nambu 
(particle-hole and spin) space and obey transport-like 
equations for excitations of energy e moving along clas- 
sical trajectories labelled by the Fermi momentum pj. 
The temperature- and phase-biased contact is, to qua- 
siclassical accuracy, in a stationary state; thus, we drop 
the dependence on t. 

The heat current through the point contact 



Ie((/>, Ti, T2) is a function of the phase difference 4> 
and the bath temperatures Ti and T2. In our ge- 
ometry, the heat current flows along the junction 
normal only. If = zl^, and is found by energy integra- 
tion and Fermi-surface averaging of the quasiclassical 
Keldysh Green's function. Current conservation allows 
us to write the heat current in terms of functions 
locally at the junction on trajectories with positive 
and negative projections of the Fermi momentum on 
the z-axis: = 6^(p/ • z > O,e;0,ri,r2) and 

^ ©^f (p^ • z < 0, e; (/), Ti, Ta), respectively. Thus, 



Ucl),Ti,T2) = ANfVf 



de 



(3) 



where Nf is the normal-state density of states at the 
Fermi level, A = tto^ is the cross-sectional area of the 
contact, and the angle brackets denote a Fermi-surface 
average, including the projection of the group velocity 
along the direction normal to the interface. 

A efficient method of computing the Green's functions 
is provided by the parametrization in terms of gener- 
alized spectral functions i<Jyj^ and distribution hmc- 
tions , which are 2x2 matrices in spin space. In 
the spin-singlet case, the spin structure of the spectral 
functions are all given by iay , while the distribution func- 
tions are proportional to the unit matrix. These scalar 
amplitudes obey Ricatti-type differential equationsi^i^SiSi 
Each Green's function can be written in terms of a set of 
Ricatti amplitudes. The retarded and advanced Green 
functions have the forms 



=F7ri / 1 



R.A 
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R.A 



-2iay^^-'^ 



R,A 



2i(JyJ 

-1 ^;yR.A^R,A 



(4) 



where (^^'^ = 1-^-jR,a~r,a^ while the the Keldysh Green 
function has the form 



(5) 



+ x'^-f^Y 

r,K~R z.K~A 



The advanced amplitudes are related to the retarded ones 
through the symmetry 7"^ = —(7'^)*. We will also make 
use of the conjugation symmetrjiiS 



q{Vf,z,e) = (7(-p/,z,-e)* 



(6) 



For each trajectory, the spectral functions and distri- 
bution functions arc found by integrating the correspond- 
ing Ricatti equations with initial condition either in the 
bulk or at the interface, depending on the stability prop- 
erties of the differential equation. We follow the nota- 
tion in Ref. 0; quantities denoted by lower case letters 
are computed by integrating the Ricatti equations with 
initial conditions in the bulk, while quantities denoted 
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with upper case letters are computed by integrating the 
Ricatti equations with initial conditions at the interface. 
According to this convention, for each of the four trajec- 
tories in Fig. n] we change lower case letters in Eqs. I4I5I 
to upper case letters; the list of amplitudes defining each 
Keldysh Green function for the scattering trajectories are 
shown in Fig. ^ 

In the bulk region, the amplitudes are given by their 
equilibrium values 



TABLE L The Andreev reflection probabilities. The denom- 
inators are listed in Table UTI 
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!,i4>j 



A,(p/)e-'^^- 
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(7) 



xf(p/,e) = (l-|7f(p/,e)|^)tanh — 



2Ti 



if(P/,£) 



(l-|7f(p/,e)ntanh 



2T, 



where the index j = 1, 2 refers to quantities in su- 
perconductor one (two). Both and the mean field 
gap Aj formally include renormalization effects (self- 
energies) from elastic and inelastic scattering. In the 
following we consider the ballistic (clean) case and de- 
fer the discussion of impurity scattering to Sec. IVII We 
also assume that the inelastic scattering rate, Fj^, is small 
compared to all relevant energy scales. We can then write 
£^ = e-|-iFi„ e-|-iO+, and Aj denotes the temperature- 
dependent weak-coupling gap amplitude. 

The lower case amplitudes are found by integrating the 
Ricatti equations from the bulk to the interface, with the 
initial conditions in Eqs.[7| The initial conditions for the 
upper case Ricatti amplitudes at the contact are then 
found from Zaitsev's non-linear boundary conditionSfS^ 
which in terms of Ricatti amplitudes are reduced to linear 
boundary conditions^ For the distribution functions we 
use the notation in Ref. [23 and obtain. 



(8) 



In terms of the following particle-hole spinor notation. 
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(a| = (l, -laya*), (9) 



the Green's functions at the junction (z = 0) can be 
written in a compact form 

0f+ = [xf |r,e)(r.e| +lf ri|7f)(7f In' 



K 



©f_ = 



K 
2+ 



2711 
2Tri 

-— [x^T,\reH){reH\r,+X^\^^){^^\] 



[iffi|r,,)(r,,|n+Xf|7f)(7f|] 

'x^\fue){fhe\+X^h\l!^){l^\h 



(10) 



where Nk = |CfeP for k = 1...4 with 



1 + l2rhe : 



Ci = l+7inie, Ca 

C2 = 1 + ll Teh , C4 = 1 + J^feh . 



(11) 



The Andreev reflection amplitudes are listed in 
Tabled and the scattering probabilities. Tap and Rap, 
arc listed in Table UTI The notation is chosen such that 
(Rap) denotes the probability of transmission (re- 
flection) of a particle of type /? to a particle of type a. 
Quantities with a bar denote transmission from right to 
left and reflection on the right side, while those without 
bar denote transmission from left to right and reflection 
on the left side. 

We use the above interface Green functions to compute 
the heat current via Eq.|31 There are four contributions, 
one for each incoming distribution function in Eq. i = 
1 — 4 for xf' ,Xi ,X2 ,X2 , corresponding to electron- like 
and hole-like quasiparticles injected from the left {i = 
1 — 2) and right (i = 3 — 4) reservoirs, 



J —00 \ .-I 



(12) 



We can express the spectral current in terms of trans- 
mission probabilities, T^, for the distribution functions. 
The most direct way of achieving this is to make use of 
the unitarity of the scattering matrix for the junction, 
which leads to the relation 



Tr 



6j(p/,z = 0-,e)-©f(p/,z = 0-,e) 
0f (p/, z = 0+,e) - ©f (p/,z = 0+, e) 



Tr 



(13) 



Thus, we compute jl and at 2; = 0^, while and 
are computed at z = 0^. The heat current spectral 
densities can then be written as 
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TABLE II: Scattering probabilities for distribution functions in the stationary SIS junction setup. 



T„„=D|l+7f72«lV!Bi|' 






Bi = 1 + i?72«72« + D^^'y^ 


Tee=7?|l+7f72^|V|B2p 




R,,^R\l+j^j^\y\B,\^ 


B2 = 1 + R-y^J^ + D-y^j^ 






Rhh=R\l + i?-^?\V\B:i\^ 


Bs = 1 + i?7f 7f + 0^2 71' 


T,, =73|l + 72«7f|V|B4p 






B4 = 1 + i?7f 7f + £>7l'7f 



i,i(p/, = 0+, e; 0, Ti, T2) - -ANfVjex^ 
jUPf,z = 0+, e; 0, Ti, T2) = -^iV/v^eif 
i3(p^, z = 0-, e; 0, Ti, T2) = -^iV/z;/ex^^ 
j,^(p/, z = 0", e; 0, Ti, T2) = -^iV/z;/e5;^ 



[-T,e(l - |7f HA^r' - - |7f nA^2"i] , 



Note that each spectral current density vanishes in the 
subgap region, i.e. for |e| < max(Ai, A2) the factor 1 — 
1 7^ 1 2 = The heat current is only carried by continuum 
energy quasiparticles, and thus, we only consider |e| > 
niax(Ai, A2) from here on. Using the symmetry in Eq.|Bl 
we obtain 

j2(p^, z = 0+, e) = jli^pf, z = 0+, -e), 
j,(p/,z = ,e) j,(-p/,z = ,-e), 

which implies that hole- like quasiparticles carry the same 
amount of heat as the electron-like quasiparticles (after 
energy integration and Fermi-surface averaging). 

We introduce transmission coefficients (script S)'s) by 
combining the scattering probabilities for x^"^ distribution 
functions (big T"s) with the spectral renormalization fac- 
tors [the factors, (1— |7''^P), and the denominators N~'^]. 
We then obtain. 



jl = -ANfVfe tanh— [Dee{Pf,e) + ^he{Pf,e)\ 
fe = +ANfVfe tanh^ [S)ee(p/, e) + ©/^^(p/, e) 



(16) 



The transmission coefficient for electron-like quasiparti- 
cles remaining electron-like, See, has the form 



2)ee(P/,e) = (1 - |7fT)ree(l - |72''f )^^' 

D{i-\i^?){i-\im^+i^i^\'' (17) 



while the transmission coefficient for electron-like quasi- 
particles with branch conversion to hole-like quasiparti- 
cles, Dhe, is 

D/ie(P/,e) = 7^ • (1») 



right to left are, 

D{l-\lf?){l~\l?f)\^ + lMf 



®ee(P/,e) 



D/ie(p/,e) 



fliJ(l-|72^n(l-|7fni7i^-72^| 



(19) 



The common denominator is given by 

Z=l + i?(7i%« + 72''7i') 

+ ^(7i''72''^ + 7i''72'') + 7i''7i''72''72''- 

We immediately see that 

£'/ie(p/,e) = D/ie(P/,e), 
while the symmetry given by Eq. |B| leads to 

See(P/,e) = Dee(-P/,-e)- 



(20) 



(21) 



(22) 



With the aid of these relations we can express the total 
heat current in a more compact form, 

/>CXD 

Ie(<^,Ti,r2) = 2 / deje(^ = O,e;0,ri,r2), (23) 

J max( A) 

where the factor 2 reflects the symmetry of the spectral 
current under —e — > e. The heat-current spectral den- 
sity is expressed in terms of the Fermi-surface averaged 
transmission coefficient, 



je (z = 0, e; 0, Ti , Ts) = -2ANfV / e ( J) (p/ , e, 0; Ti , Ts) 

(24) 



tanh — - — tanh — — 
2ri 2r2 



The corresponding coefficients for transmission from 



2)(p/,e,0) =l)ee(p/,e,0)+£>he(p/,e,0), (25) 
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where the factor two is due to the electron-hole symmetry 
[c.f. Eq. Equation [^satisfies the symmetries 



i,(-0,ri,T2) = i,(0,ri,T2), 



(26) 



and Eas. l23l2H arc the results for non- linear heat-current 
response to a temperature and phase bias. These results 
hold for singlet superconductors with any orbital symme- 
try, in particular for either s-wave or d-wave symmetries. 



A. S-wave symmetry: A(p/) = A 

In the remaining sections we focus on the low-T^ su- 
perconductors, with a momentum independent (s-wave) 
order parameter. However, the order parameters for the 
two superconductors have different phases, (jjj, and may 
also have different order parameter amplitudes, e.g. as 
a result of a finite temperature bias, Aj ~ A{Tj). The 
retarded Ricatti amplitudes in the bulk are then momen- 
tum independent. Moreover, since the s-wave order pa- 
rameter is constant in space for the point contact, the 
lower case Ricatti amplitudes are given by their bulk val- 
ues along the whole trajectory. We then obtain analytic 
expressions for the effective transmission coefficients 



®ee(P/,e) 



S/!.e(P/,e) = 



2j^66(£' + 66- A1A2 cos(/)) 
[De^ + (1 + i?)ae2 - DAi A2 cos ( 

2i^jja6(£'-ei6- A1A2 cosq 
[De2 + (1 + _ Z?Ai A2 cos ( 



(27) 



where = sgn(e)4/ e 



B. Linear response 

In linear response we write Ti ~ T and T2 — T + 5T . 
Then to lowest order in 5T we obtain. 



h{<P,T)=~K{^,T) 5T, 



where the heat conductance has the form 



(28) 



df 



K{cj>,T) = AAj^ &AA(e)K(e)] (!)(£,</., T)) j , 

(29) 

with the transmission coefficient, 

S)(p/, e, <j>, T) = S)ee(p/, e, 0, T) + S)„e(p/, e, 0, T) , (30) 



See(P/,e,0,T) 



(e2 - A2)(e2 



£2- A2(l-Dsin2|) 
(e2 - A2) A2 sin^ | 



(31) 



- A2(l-Dsin2|) 



Equation Hi is an intuitive form of the heat conduc- 
tance that resembles the bulk thermal conductivity. The 
conductance is expressed in terms of the bulk quasipar- 
ticle density of states, A^(e) = A^/|e|/\/e2"^^"A2, and the 
energy current carried by these quasiparticles, [evg{e)], 
where Vg{e) — Vj\Je^ — A2/|e| is the group velocity of a 



bulk excitation. Note that the product, Mvg 



N 



is energy independent. The backscattering at the junc- 
tion (resulting from the Sharvin resistance and the trans- 
parency D < I), which limits the heat conductance, cor- 
responds to the elastic mean free path due to impurities 
in the bulk. 



The linear response results in Eqs. 128131)1 for the 

heat conductance were reported in Ref. ,4, for a p/- 
independent normal-state transmission probability, D. 
In Sec. IIVI we extend the analysis of the linear re- 
sponse limit and examine the non-linear thermal re- 
sponse, Eqs.ESd in Sec.|V| 



IV. THERMAL CONDUCTANCE 

The phase modulation of the heat conductance, 
k((^, T), is determined by the transmission coefficient, 
D{e,4>,T). The energy dependence of the transmission 
coefficient reflects features of the phase-dependent local 
density of states (LDOS) of the junction, J\fj{e,4>,T). 
The LDOS at the contact is defined as A/j(p/, e, 0, T) = 
-^Im(5fi(p/,z = 0",e, 0,T), where (Sf^ is the diag- 
onal (11) component of the retarded Green's function, 
which has the following form at the junction. 



6fi(p/,z = 0-,e) _ eVA2 - e2 - |sgn(p/ ■ z)A'^Dsm, 



e2- A2(l-Dsin^f) 



(32) 



The resulting LDOS has both a bound-state (|e| < A) 
and continuum (|e| > A) spectrum given by 
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AO(p^,e,(/.,r) 



e(|e|-A) 



AVI 



+ 7r(5(e - es) 



e2- A2(l-Dsin2(0/2)) 

'e|VA2 — — ^sgii(p/ • z)A2£)siii( 



2|e| 



7r(5(e + 



|e|VA2 - £2 + 5Sgn(p/ • z)A2Dsin(; 

M 





(d) 



it)=0 
(t)=0.47i 
it)=0.6Jt 
(t)=n: 



FIG. 2; The LDOS at the junction, Nj{e,(f)), and the trans- 
mission coefficient, S(e, (^), at fixed momentum p/ with D = 
D{pf). (a)-(b) DifTerent transparencies at pliase difference 
(j) = n. (c)-(d) Different phase differences at transparency 
D = 0.2. Nj is suppressed by the formation of the ABS, 
as is tuned from to tt. For high transparency D is sup- 
pressed, but for low transparency S) contains a resonance at 
€r — A + iADsin^ ^, at which Z){er,(f>) = 1/2. In all cases. 



T — 0.72Tc, and for clarity we added a small width, ~ 10 
to the bound-states. 



The spectrum is shown m Fig. [21 as a function of the 
phase bias ^) and as a function of barrier transparency 
for (j) = n ^^). Of particular importance is the forma- 
tion of a pair of Andreev bound states (ABS) with phase 
dispersion 



±Aa/1 - Dsin^ 



2 ' 



(34) 



and the impact of the ABS on the continuum spectrum. 
In addition, we note the following characteristics of the 
LDOS: 

1. In the absence of phase bias, = 0, the LDOS 
reduces to the BCS density of states in the bulk. 

2. Under a phase bias, (p ^ 0, Andreev bound states 
are formed, with spectral weight drawn from the 
continuum near e = A. 

3. The ABS are weakly bound and close to the con- 
tinuum edge for low transparency. For high trans- 
parency these states are more strongly bound, and 
may lie well below the gap. 




FIG. 3; The thermal conductance as a function of <j) and bar- 
rier transparency D. The thermal conductance is normalized 
for each D by its value at = 0. 



4. For (f> = TT the ABS is the closest to the Fermi 
level (at es ~ A\/i?), and consequently the spec- 
tral weight of continuum excitations is reduced the 
most at this phase bias. 

These characteristics of the LDOS have direct conse- 
quences for the properties of the transmission coefhcient, 
D, and thus for the heat conductance. For = the 
transmission coefficient becomes independent of energy 
and reduces to the normal-state transmission probabil- 
ity: D{e,0) = D. Thus, the thermal conductance of the 
point contact at = 0, 



k(0 = 0) 



NfVf {D)A 

rp2 



dee^scch^ 



e 

2r ' 



(35) 



is reduced compared to the normal-state conductance 
(Eq.OHwith A = 0), AtN = ^NfVf {D)AT, by opening 
the superconducting gap. Note that the temperature de- 
pendence of the ratio, k(0 = 0)/kn('7c), is equivalent to 
that of the normalized bulk thermal conductivity of an s- 
wave BCS superconductor. For 7^ 0, the transmission 
coefficient is strongly dependent on energy and phase, 
which is the source of the phase-modulation of the ther- 
mal conductance shown in Fig. |21 For high transparency 
{D > 0.7), the suppression of the continuum density of 
states associated with the formation of the ABS well be- 
low the gap is reflected in a suppressed transmission co- 
efficient and a reduced heat conductance compared to 
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= 0. Thus, for high transparency as (j) is tuned from 
to TT, the conductance is suppressed as shown in Fig. O 
In the hmit D ^ 1 only direct transmission is possible, 
Dee(e» = (e^ - A^)/(e^ - A^cos^^)); branch conver- 
sion processes vanish, Sl/ie = 0, and we recover the heat 
conductance of a pinhole junction 



(36) 



obtained by Kulik and OmelyanchukiSi 

In the case of low transparency, D ^ 1, and 7^ 0, the 
Andreev bound states lie just below the gap edge. The 
transmission coefficient, S)(e, 0), exhibits a resonance just 
above the gap edge, as shown in Fig. [5| This resonance 
is a consequence of the shallow ABS. We sec this fact 
clearly in Fig. |5Jc)-(d): the continuum LDOS is reduced 
as (j) is changed from to tt, while the resonance in the 
transmission coefficient is enhanced (the area under the 
curve is enhanced). The resonance energy is obtained by 
writing = A+(5e,., where 5er ^ A, in the expression for 
the transmission coefficient. Eg. 1301 From the extremum 
condition, dTl/dSer = 0, we obtain 



A ( 1 + -Dsin^ - 
2 2 




(37) 



FIG. 4; The temperature dependence of phase-modulation 
of the thermal conductance for (j) = n, normalized by the 
conductance at Tc, k,n = AN fV fTc (D) /3. Shown for 
comparison is the normalized conductance for = 0, which 
is independent of transparency. 



is a consequence resonant transmission of heat carrying 
quasiparticles, and is characteristic of intermediate- to 
low-transparency Joscphson point-contact junctions. 



At the resonance energy the transmission coefficient is 
independent of transparency and takes the value 



S)(er 



(38) 



The resonant transmission of quasiparticles with ener- 
gies e « er leads to an increase in the thermal conduc- 
tance as the phase bias is tuned from to tt, as shown 
in Fig. 121 Observation of the phase modulation may be 
accomplished with a point-contact Josephson junction in 
a SQUID geometry. In this case the phase is tuned by 
the flux, $, threading the SQUID, = 27r($/$o), where 
$0 = hc/2e is the superconducting flux quantum. 

Another important consequence of the resonance is the 
increase of the thermal conductance for (f> — tt compared 
to the normal-state conductance at Tc, as the temper- 
ature is lowered below Tc, see Fig. 0] This peak effect 
is large and persists over a broad temperature range, 
0.5Tc ^ T < Tc, for moderate to low transparencies, 
D ~ 0.2 to £> < 1. 

For very low temperatures, T ^ A, the sharpening 
of the distribution functions at the Fermi level, com- 
bined with the gap in the continuum spectrum leads 
to an exponentially small thermal conductance. The 
reduction of the low-temperature heat conductance for 
normal-superconducting (NS) proximity structures is 
well known from Andreev's work on heat reflection at 
NS interfaces^ However, the increase in the thermal 
conductance compared to the normal state in the in- 
termediate to high temperature range is not related to 
Andreev reflection at high-transparency interfaces, but 



A. The tunnelling limit: D ^ 

Here we examine in more detail the phase dependence 
of the thermal conductance in Eq. 1291 in the tunnelling 
limit, _D <C 1. We neglect the momentum dependence 
and set D{pf) = D. For D 1, when we expand 
the transmission coefficient, S)(e, 0), to leading order 
in D there is prefactor proportional to D, and a non- 
perturbative dependence on D that is responsible for the 
resonance peak at ej. The resonance leads to a non- 
analytic contribution to the conductance proportional 
_D ln_D that dominates in the limit D 0. Consequently 
up to order D the thermal conductance k(0) contains a 
term proportional to sin^ ^ , as well as a nonanalytic term 
proportional to sin^ ^ ln(sin^ ^), 

k(0) = ko - ki sin^ ^ ln(sin^ ^) + K2 sin^ ^ 

where kq is given by Eq. 1351 and 

NfVfDAA\ ,A 
^1 - sech 



(39) 



K2 



NfVfDAA^ 
4r2 



4T , A , 

— — ( 1 — tanh — ) -|- c 
A^ 2T' 



{1 + IiiD)ki, 



dxx In X 



sinh 



2T/A 



sech 



2^/^+1 
2T/A 



(x2 + l)cosh3^ (a;2 + 1)3/2 
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D=0.01 
D=0.005 



T=0.5 




D=0.01 
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K, term 
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continuum 
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1 



FIG. 6: Three regions in energy: sub-gap energies below the 
gap of both superconducting leads, semi-continuum energies 
between the two gap energies, and the true continuum. Only 
the true continuum states transport heat. 



FIG. 5: The ratio ki/k2 as a function of temperature. The 
inset shows the phase dependence of the thermal conductance, 
the contributions from terms proportional to ki and K2 are 
plotted separately. 



The phase modulation of the thermal conductance is 
shown in the inset of Fig. |31 the term K2 sin^ ^ and 



the term — ki sin -I In (sin ^) are plotted separately for 
T = 0.5Tc and D = 0.01. The relative importance of 
the two terms is shown for D = 0.01 and D = 0.005 as 
a function of temperature. The ratio ki/k2 increases as 
temperature is lowered. 

We compare Eq. 123 with the result obtaine d with the 
tunnelling Hamiltonian (tH) method (see Refs. 1 121131 and 
Appendix). According to the tH calculation, the heat 
transport through a tunnel junctions is given by Ea. lA.3l 
in the Appendix. In the linear response limit, ST —^ 0, 
Ai = A2, and Ea. lA. 31 reduces to 



87riV||Tp 



de e- 



cos ( 



dT 



(40) 



The integral in Eq.^Ujis divergent due to the singularity 
at e = A. The divergence is unphysical and indicative 
of the failure of low-order perturbation theory within the 
tH method for the linear response. What is missing in 
Eq. 1401 is the correction to the spectrum - which enters 
the denominator - that results from the small, but finite, 
barrier transparency. This non-perturbative correction 
is fully taken into account in the quasiclassical Green 
function method. The singularity of the DOS at e = A 
is removed by the formation of the ABS, and the final 
result, Eqs. I29I30I contain no divergence in the linear 
response limit. 

In earlier treatments the divergence entering the tH re- 
sult was regulated by introducing an ad-hoc cutoffjiSii^ 
for example by requiring that the two gaps have differ- 
ent magnitudes, Ai 7^ A2. Under these circumstances, 
Eq. IA.3I predicts the thermal conductance has the form 
k(0) = Kq — ^2 cos (f), with both k'q and K'2 proportional to 
|Tp oc D. However, as shown in Eq. I^Hlthe thermal con- 
ductance of a tunnel junction contains non-perturbative 



corrections to the cosine dependence, i.e. the leading 
order correction is —ki sin^ j ln(sin^ ^). Also the mag- 
nitude of K2 contains a term proportional to DlnD. 



V. NON-LINEAR RESPONSE 

When we consider the non-linear response one impor- 
tant effect is that the order parameters on the two sides 
have different magnitudes. 



Ai = A(r), 

A2 - A(r + AT) < Ai. 



(41) 



The density of states is modified accordingly. There are 
three important regions of the local excitation spectrum 
(c.f. Fig.EJ: 

1. |e| < A2, sub-gap spectrum with bound states, 

2. A2 < |e| < Ai, semi-continuum spectrum, 

3. Ai < |e|, true continuum spectrum. 

Only the states in the true continuum carry heat. In 
Fig. [7| we plot the local density of states at the junc- 
tion for three different phase differences. The differ- 
ence between the gap magnitudes becomes large as AT 
is increased, the bound states move closer to the semi- 
continuum edge. At the same time, the true continuum 
edge is further separated from the ABS. Consequently, 
the true continuum is less affected by the bound states, 
and therefore the phase bias. 

The Andreev bound states also have an impact on the 
heat current in the non-linear response, via a mirror ef- 
fect in the continuum which introduces a resonance in 
the transmission coefficient. However, in the non-linear 
response, the true continuum is shielded from the bound 
state by the semi-continuum. This weakens the reso- 
nance and suppresses the signatures of resonant trans- 
mission in the heat current. The transmission coefficient 
is plotted in Fig. |S{c). For low transparency, the reso- 
nance is weakened when AT is increased. More specifi- 
cally, the area under the transmission curve, which is the 
relevant quantity for the heat current, is reduced. For 
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N(E. ^=71/2) 



E/A, 



0.6 0.7 o.« o.y I 

e/A, 



FIG. 7: Local density of states at the junction sd, z — 
in the non-linear response, for transparency D = 0.5 and 
temperature T = 0.72Tc. 



high transparency, the opposite happens: the transmis- 
sion coefficients is less suppressed near the gap edge as 
AT is increased (not shown). 

In addition to the spectral change found in D, the 
change of the thermal occupation factors for the two 
reservoirs influences the heat transfer. We have checked 
that the main effect on the heat current is due to the 
spectral change, but the change in the occupation fac- 
tors reduces the overall effect of a phase bias on the heat 
current. The total heat transfer is larger for a larger 
temperature bias, but we eliminate this scale factor by 
normalizing the heat current by the corresponding heat 
current at zero phase bias, or at Tc- 

The phase-modulation of the heat current is plotted 
in Fig.|SJb), while the temperature dependence is shown 
in Fig. IHl^a). Since the angle dependence of the barrier 
transparency does not qualitatively affect our results, we 
use the angle-independent model. The phase-modulation 
is weakened, and the peak effect near Tc for low trans- 
parency and phase difference (f) = n is reduced. However, 
the non-linearity does not remove or drastically alter the 
resonance effects, and the main features found in linear 
response persist for finite temperature bias. 
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FIG. 8; (a) Temperature dependence and (b) phase depen- 
dence of the heat current in the non-linear response, (c) 
Transmission coefficient in the non-linear response at phase 
bias <j} — n. In (b)-(c) the temperature is T = 0.72Tc, and in 
all cases the transparency of the junction is _D = 0.1. 



by 
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J deeTvigVg)''. 



(43) 



To calculate the heat current through the point contact 
with dirty leads, we need the boundary conditions for the 
isotropic propagator, g, at an interface. Such boundary 
conditions were derived by Nazarov in the formulation 
of the circuit theory of diffusive hetero-structuresiS& In 
this model the contact is described by a set of trans- 
mission eigenvalues, {-D„}, or cquivalently a distribution 
function, p{D). The boundary conditions are 



^ANf2Vf2hg2dzg2 



dDp{D)i, 



(44) 



(45) 



VI. DIFFUSIVE CASE 

In this section we discuss the case where the two leads 
are dirty superconductors, i.e. the elastic mean free 
path, ^ ^ ^A- In dirty superconductors, ©(p/) is nearly 
isotropic, and 



g(R;e) = — 

—in 



Pf 



©(p/,R;e,) 



47r yo 

'(42) 

satisfies the Usadel diffusion-type equationi2^ The 
isotropic propagator determines most of the physical ob- 
scrvablcs. For example, the heat current density is given 



where the Keldysh matrix J is defined as 

2D[g2,gi] 
4 + D{{g2,gi}-'2) ' 



(46) 



Nfi is the normal-state density of states at the Fermi level 
of lead i ~ 1,2, Vfi is the Fermi velocity of lead i, and A is 
the cross-sectional area of the contact interface. Eq. ^| 
is the conservation of spectral current, and Eqs. I43I45I 
imply that the total heat current through the contact 
can be calculated from /, 



dDp{D) Tr(/^) . 



(47) 
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The Kcldysh component of / can be worked out to heS^ 

2~D + D{gf,g^}/2 2-D + D{gf,g^}/2 
2(2- D + D{g^, g^}/2) {2~D + D{gt,g^}l2) 

where (1 ^ 2) means the exchange of index 1 and 2 in 
all three terms. 

For point contacts, we approximate g^'^'^ with their 
bulk values, e.g. gf = {gf - gf )tanh{e/2T,). We fo- 
cus on the linear response, assume that superconducting 
leads 1 and 2 have the same gap and evaluate the con- 
ductance in the limit ST ^ 0. After some algebra we 
find 

= '2^ I '^^P^^^i'" - 

£)(e2-A2cos|) + i^i?A2sin2(,/>/2) , , 
X — ^ ^ . (48) 

Notice that with the replacement, 

J dDp{D)[...) => 2TrANfVf{...) , (49) 

the result for the conductance with diffusive leads, Eg. 1481 
in the limit of a narrow distribution of transmission bar- 
riers, is the exactly the same as Eqs. I28I30I for the con- 
ductance obtained in the ballistic limit. 

The near equivalence of the conductances for the dif- 
fusive and ballistic limits is related in part to the small 
expansion parameter, o/^a- In this limit pairbreaking ef- 
fects resulting from the current flow in the vicinity of the 
junction can be neglected, and as a result the impurity 
renormalization of the excitation energy and order pa- 
rameter cancel for s-wave superconductors. The resulting 
thermal conductance of point contact is then insensitive 
to impurity scattering in the leads. This result is not 
expected to hold for larger area contacts with a > 
These junctions are beyond the scope of this work. 

VII. CONCLUSIONS 

We have derived a general expression for the phase- 
and temperature-dependent thermal current through 
small Josephson weak links. The results are valid for 
arbitrary transparency of the junction. In the ballis- 
tic limit, we obtained an expression for the heat cur- 
rent in terms of a transmission coefficient, 2)(e, 4>,Ti,T2), 
for heat transport by continuum quasiparticle states. 
The transmission coefficient includes direct transmission 
processes, Dee, and transmission processes with branch 
conversion, Dhe- The phase modulation, temperature 
dependence, and the dependence on the barrier trans- 
parency of the heat current can be understood intuitively 
in terms of the properties of the transmission coefhcient. 



In linear response, for high transparencies, the suppres- 
sion of the continuum density of states associated with 
the formation of a low-energy Andreev bound state leads 
to a suppression of the transmission coefficient, and the 
thermal conductance, as the phase difference across the 
junction is tuned from to tt. As a consequence, the ther- 
mal conductance drops below the conductance for ^ 
for temperatures below Tc- 

However, for intermediate to low transparencies, a res- 
onance develops in the transmission coefficient as the 
phase difference is tuned from to tt. Consequently, the 
conductance is larger at (/) = n than at (/) = 0. The 
transmission resonance is due to the ABS being close to 
the gap edge. This effect is similar to transmission reso- 
nances found in wave mechanics for a quantum well with 
a shallow bound state just below the continuum edge. 
For (f> = TT the resonance in the heat current leads to 
an increase in the thermal conductance compared to the 
normal-state conductance at Tc, over a broad tempera- 
ture region below Tc. 

In the low-transparency limit, Z? ^ 1, we derived an 
analytic expression for the heat conductance and found 
that there are non-analytic corrections, of type D\nD, 
to the usual linear in D term. Also the phase dependence 
contains non-analytic terms of the form sin^ ^ ln(sin^ ^). 
The first-order tunnel Hamiltonian method gives a di- 
vergent result in linear response, which is regulated by 
the formation of an Andreev bound state from the con- 
tinuum states near e = A. A similar divergence occurs 
in the sub-harmonic gap structure (SGS) of the current- 
voltage characteristics of a superconducting tunnel junc- 
tion when the sub-gap structure is computed to finite 
order in perturbation theory within the tunnel Hamilto- 
nian method^ In the case of the SGS, summation to 
infinite order within a wave-function method)22iSS or the 
tH method,^" regulates the divergences. Such a summa- 
tion is possible for the thermal conductance within the 
tH method. However, the quasiclassical Green's func- 
tion method with the interface boundary conditions in- 
cludes the bound-state spectrum so unphysical diver- 
gences never appear. 

We also studied the case when the superconducting 
leads are in the diffusive limit. Based on the junction 
boundary conditions developed by Nazarov, we found 
that the heat conductance has the same form as in the 
ballistic case. 

Finally, in the non-linear response we find that the res- 
onance in the transmission coefficient is reduced by the 
presence of the semi-continuum spectrum, which shields 
the extended continuum energy quasiparticle states car- 
rying heat from the ABS. The effects found in the linear 
response are reduced, but not dramatically. Thus, we 
conclude that it is not essential to be in the linear re- 
sponse limit in order to observe the resonance effects in 
the heat current. 
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APPENDIX: PERTURBATION THEORY 

We review the first-order perturbation calculation of 
the heat current for low-transparency junctions based on 
the tunnelling Hamiltonian (tH) methodi^ Consider two 
superconductors, labelled by L and R, weakly coupled by 
an insulating layer. They are assumed to be described by 
the Hamiltonian 

H total = Hq + Ht, 

where Hq is the sum of the BCS reduced Hamiltonian of 
superconductors L and R, 

Ho = Hl + -ff_R, 



p,(7 p,p' ,a 



k.(T k.k' .(7 

and Ht describes the tunnelling processes, 
Ht=Y [Tkpa\f. ,^Cpa + h.c]. 

k,p,(T 

The two superconductors are also assumed to be at the 
same chemical potential = ^J■R^ but at different tem- 
peratures Tl ^Tr. 

The momentum index, p, and the annihilation opera- 
tor, Cp.cr, are reserved for superconductor L, while index, 
k, and operator, ak,a, are reserved for superconductor 
R. The spin state is labelled by a, and Tkp is the tun- 
nelling matrix element. With the BCS approximation, 
Vppi = A < 0, the order parameter is defined as 



The heat current operator is defined as 



I. = — 



dt 



= -i[Htotah Ql] = -i[HT-, Ql], 



Ql = Hl - fJ-Ly^ 4 



It is straightforward to work out the commutator 
[Ht, Ql]- Within the mean field approximation, the pair 
annihilation operator c-q^-^Cq^a reduces to its average 
value, and we find 



t 



kp^Lak,a^-p,~ 



h.c. 



(A.l) 



with ^p = Sp- fiL- 

To calculate the ensemble average of the operator I^, 
we use the interaction picture and treat Ht as a per- 
turbation to Hq. According to first order perturbation 
theory, the heat current has the form 



where Ie{t) and TiT{t') are the operators in the inter- 
action picture corresponding to /<; and Ht- respectively. 
The ensemble average (...)q is defined by iJp, and t] = 0"*". 
The commutator [Xe(t), 7iT(t')] contains various correla- 
tion functions, 

/oo 

where we define Cij(t,t') = ([0,(t), T'j (t')])o with 
Oi(t) = Tkp^P^iAt)cpAt). 

k.p.a 

= -Y TkpALal^it)clj,^_,{t), 
^iW = E Tkpaijt)cp,,it), 

V2{t)=Vl{t), 

and ak^a{t) and Cp^a{t) are annihilation operators in the 
interaction picture. We introduce the following single- 
particle correlation functions, 

G<,{t,t') ^ M{cl,{t')c^At))o, 
F>Jt,t') = -i{c^p,_,{t)cp,,{t'))^, 

F^A^,t') = +i {cpAt')c-p.-At))o, 
F<At,t')^+^(cl^,_AtVp.At)] 



/ J '^p,<T'^P,Cr- 

p,a 



The correlation functions, , can then be decomposed 
into G's and F's with the aid of Wick's theorem. For 
example, 

Cii[t,t') = - E TkpT^k-pS,p 

k,p,a 

X [FlAt^ t')F<^{t\ t) - F<^[t, t')F>Jt\ t)]. 
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The contribution of Cu to the heat current is 

/oo 
di'e''*'e(i-t')Cii(t,t') 
-oo 

Ef duj f duj' 

k,p 

CJ' — CiJ — IT] 

The factor 2 conies from the sum over spin, and A and 
B are spectral functions, with 

G^{^) = Tj/L(Tw)^b,w), 

if H = t«/l(Tc^)B(p,w). 

Here we introduced the Fermi function jhij^^ — {&^^^^ + 
1)^^, and the spectral functions 



A(j>M = ^[(1 + ^)b{Lo - Ep) + (1 - ^)5{lo + Ep)], 

£jp typ 

B{p, Lo) = TT^iSiu - Ep) - 5[u + Ep% 
tjp 

with Ep = + A|^)^/^. Analogous calculations for the 
other correlation functions follow. Summing up all the 
contributions for the heat current gives 

k.p 

f ,2^ A{p,^)A{k,u') ^ ,^ 2a B{-p,u;)Aik,u;') 

i Mfcpl -. — h \J-kp\ Al- 



W' — + 177 

LU' — LU — IT] 

A{~p,uj)B{k,u') - 

+-l-kpJ--k-p^L ; : 

Uj' — UJ — ITj 



Ul' — OJ + IT] 



(A.2) 



In Eq. IA.2I we neglect the momentum dependence of 
the tunnelling matrix elements, put T^p — > T, and change 
the momentum sum into integrals over excitation ener- 
gies Ep and Ek- The integrals over to and uj' collapse 
because of the S functions in the spectral functions A, 
B and B. Taking the imaginary part of the integrand 
yields new S functions S{Ek zk Ep), which lead to a col- 
lapse of one of the energy integrals. Finally, the heat 
current takes the form 



r + oc 



- lAr IIArI C0S( 



v/e2-|Ai|Ve'-|Ai^P 



(A.3) 



Here the phase difference, (/) ~ (pL — 4'Ri ^max = 
max[|Ai|, |Afl|], and N^^^ are the density of states at 
the Fermi level for L and R sides, respectively. 

Finally, we compare with previous calculations of the 
heat current in tunnel junctions. In Ref. ^3 only the 
kinetic energy term was included in deriving the current 
operator; the term, TfepA^aj^ n^^ -v.-a ^ Eind its Hermitian 
conjugate were missing from Ea. lA.ll In the main result 
of Guttman et al. (Eq. 9 of Ref^U.) the sign of the cos 
term is incorrect. 
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